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Motivated by recent experiments on isovalent-doped 1111 iron-based superconductors 
LaFeAsi_a;Pa;Oi_j,Fj, and the theoretical study that followed, we investigate, within the five orbital 
model, the correlation between the spin fluctuation and the superconducting transition tempera¬ 
ture, which exhibits a double dome feature upon varying the Fe-As-Fe bond angle. Around the first 
dome with higher Tc, the low energy spin fluctuation and Tc are not tightly correlated because the 
finite energy spin fluctuation also contributes to superconductivity. On the other hand, the strength 
of the low-energy spin fluctuation originating from the d^^zjyz orbital is correlated with Tc in the 
second dome with lower Tc- These calculation results are consistent with recent NMR study, and 
hence strongly suggest that the pairing in the iron-based superconductors is predominantly caused 
by the multi-orbital spin fluctuation. 

PACS numbers: 74.20.-z, 74.70.Xa, 74.25. nj 


I. INTRODUCTION 

It is now well known that the superconducting transi¬ 
tion temperature (Tc) of the iron-based superconductors 
can vary, in some cases largely, by substituting some el¬ 
ements with others. In particular cases, the Tc variation 
upon element substitution has been found to be non¬ 
monotonic. These observations have sometimes lead to a 
speculation that the Cooper pairing in the iron-based su¬ 
perconductors may involve multiple pairing mechanisms. 
One of the recent important example is the hydrogen- 
doped 1111 systemi. Namely, large amount of electrons 
can be doped by substituting O with H in LnFeAsO with 
{Ln=Gd, Sm, Ce, La). Surprisingly, it has been shown 
that the superconductivity appears even up to 40% of 
electron doping, and particularly in LaFeAs(0,H)L and 
SmFe(As,P)(0,H)^ the phase diagram exhibits a double¬ 
dome structure against the electron doping. Some of the 
present authors of the present paper have clarified the ori¬ 
gin of this double dome phase diagrarti^. Namely, while 
the first Tc dome is mainly due to the Fermi surface nest¬ 
ing, the second dome originates from the spin fluctuation 
enhanced by a peculiar real space motion of electrons 
within the d^y orbital, despite the degraded nesting. 

Another example was recently found in an isovalent 
doping 1111 system, where As is partially substituted 
with P in L7rFeAs(0,F) (Ln=Nd, Ce, La). It is known 
that increasing the phosphorous content enlarges the Fe- 
Tn-Fe (Tra=pnictogen) bond angle (or reduces the pnic- 
togen height). Therefore, it was expected, from the em¬ 
pirical Tc trend found by Lee et al^, that Tc monotoni- 
cally decreases by increasing the phosphorous content^^— . 
However, it has recently been revealed that another local 
maximum of Tc (we will call this the second Tc dome) ex¬ 
ists in the high phosphorous content regime^i^d^. More¬ 
over, for the non-fluorine-doped system, the antiferro¬ 
magnetic phase sandwiched by the two superconducting 
phases were found in the finite P doping regional — . In 


order to investigate this problem theoretically, some of 
the authors of the present paper have recently studied 
the correlation between the Fe-Tn-Fe bond angle and the 
spin-fluctuation-mediated superconductivityi^. There, it 
was shown that the eigenvalue of the linearized Eliash- 
berg’s equation. A, which can be considered as a measure 
for the superconducting Tc, and also the Stoner factor of 
antiferromagnetism, vary non-monotonically as functions 
of the bond angle, indeed exhibiting a double dome struc¬ 
ture consistent with the experiments. They concluded 
that in the small bond angle region (small P concentra¬ 
tion regime), the spin fluctuations originating from the 
dxy and d^z/yz orbitals both contribute to the supercon¬ 
ductivity, while in the large bond angle region (large P 
concentration), the dxy orbital has small contribution to 
superconductivity, and the nesting of the dxz/yz portion 
of the Fermi surface is the main origin of the supercon¬ 
ductivity. 

Experimentally, the magnitude of the low energy spin 
fluctuation is often probed by the NMR experiment. 
Namely, 1/TiT measured in the NMR experiment is pro¬ 
portional to the slope of the imaginary part of the dy¬ 
namical spin susceptibility taken in the a; ^ 0 limit. 
Therefore, the correlation between the superconducting 
Tc and the development of 1/TiT upon lowering the tem¬ 
perature in the normal state has often been discussed in 
the context of determining whether the Cooper pairing is 
mediated by the spin fluctuation or not. In fact, in some 
cases, the superconducting Tc and the development of 
1/TiT at low temperature is found to be correlatedi^^— , 
suggesting the importance of the spin fluctuation played 
in the pairing mechanism. 

As for LaFe(As,P)(0,F) mentioned above, the corre¬ 
lation between Tc and 1/TiT has been investigated in 
ref. 3 as shown in Fig[TJ It was found that Tc at finite 
phosphorous content is correlated with the development 
of 1/TiT upon lowering the temperature, indicating that 
the low energy spin fluctuation is responsible for the non- 
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FIG. 1. The NMR experiment result of 1/T\T against T for 
LaFeAsi_j;PxOi_j,Fj, and YjLai-jFeAsOy. Reprinted from 
ref. . 


monotonic Tc dependence with phosphorous doping. On 
the other hand, it was also shown in the same paper 
that in a 1111 material Lao.osYo.gsFeAsOi-y, in which 
Tc reaches a very high Tc of ^ 50K, the development 
of 1/TiT tends to be saturated in the low temperature 
regime just above Tc (FiglT)). Comparing this result to 
that for LaFe(As,P)(0,F) suggests that the Tj. and the 
development of the low energy spin fluctuation is not 
necessarily correlated, indicating that factors other than 
the low energy spin fluctuation is playing a role in the 
enhancement of the superconductivity. In fact, the ob¬ 
servation that Tc is not necessarily correlated with 1/TiT 
has already been recognized shortly after the discovery 
of the iron-based superconductors. Namely, in ref.[l^, it 
was pointed out that Tc is not largely affected by dop¬ 
ing electrons into LaFeAs(0,F), while the development 
of 1/TiT at low temperatures is suppressed prominently. 
Therefore, it seems that 1/TiT is sometimes correlated 
with Tc, while in other cases not. It is interesting to see 
whether these experimental observations can be under¬ 
stood within the framework of the spin-fluctuation medi¬ 
ated superconductivity consistently. The doping depen¬ 
dence of 1/TiT has indeed been investigated by Ikeda in 
the early daysi^d^, but the non-monotonic Tc behavior 
in LaFe(As,P)(0,F) was not known at that time. 

Given the above mentioned background, in the present 
study, we theoretically study the correlation between the 
Fe-Tn-Fe bond angle, the strength of the w 0 or finite 
energy spin fluctuation, and the superconductivity me¬ 
diated by spin fluctuation. We conclude that the exist¬ 
ing experimental results can indeed be understood within 
this framework, where the key point is the contribution 


to superconductivity coming from finite energy spin fluc¬ 
tuation that has small contribution to NMR I/TiT. 


II. FORMULATION 
A. Model construction 

We first perform first principles band calculation of 
LaFeAsO using the WIEN2k packagei^. We adopt hypo¬ 
thetical lattice structures of LaFeAsO, where the Fe-As- 
Fe bond angle a is varied, while fixing the Fe-As bond 
length. From the first principles band calculation, using 
the maximally localized Wannier orbitals22r22 and fol¬ 
lowing ref. [1^ , we construct the five orbital tight-binding 
model in the unfolded Brillouin zone in the following form 

i^j <7 

, where c| creates an electron with spin a on the /i, v 
th orbital at site i,j, and are the hopping integrals. 
We define the band filling n as the number of electrons 
per number of sites (e.g., n = 10 for full hlling). On top 
of the non-interacting part iLg, we take into account the 
electron-electron interaction in the following form : 

Hi - EE EE^^m. 

i fi>L' era' 

~ E/ + E/ 

where the standard interaction terms that comprise the 
intra-orbital Coulomb interaction U, the inter-orbital 
Coulomb interaction U', the Hund’s coupling J, and the 
pair-hopping J' are considered. The magnitude of the 
interactions are taken to be dependent on the orbital. 


B. FLuctuation EXchange approximation (FLEX) 


In the present study, we adopt The fluctuation ex¬ 
change (FLEX) approximation^iiS^ to the multi-orbital 
Hubbard model. The renormalized Green’s function G is 
given as 

G = (Gq-i - E)-i 


where Gg is the Green’s function of the non-interacting 
electrons and E is the self energy. In FLEX, bubble and 
ladder type diagrams are considered in the self energy 
calculation, which is determined self-consistently in the 
Baym-Kadanoff sense in order to have the conservation 
laws satisfied. The spin and charge susceptibilities ma¬ 
trices Xs{<l) and Xc{<l)^ respectively, are described in this 
approximation by the following equations. 


xM 

Xc{q) 


xo{q) 

1 - Sxo{q) 

xo(q) 

i + C'Ao(g) 
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where S and C are orbital dependent interaction vertex 
matrices, q = (q, where q is the wave vector and 
iojn = (2n + l)'i:kBT is the Matsubara frequency. The 
irreducible susceptibility matrix elements is 

given using the renormalized Green’s function as 

k 

where U (i = 1,2, ■ ■ ■ ,5) are orbital indices. 

After Xsiq) is obtained on the imaginary frequency 
axis, we perform analytic continuation exploiting the 
Fade approximation to obtain the spin susceptibility on 
the real frequency axis. We define quantities that mea¬ 
sure the strength of the intra orbital spin fluctuation as 


k 





Here, the intraorbital spin susceptibility of the /i-th or¬ 
bital is Xs = r''(a;) measures the strength of the 

spin fluctuation within the /i-th orbital at a certain fre¬ 
quency uj, while 7 ^ is a measure of the intraorbital spin 
fluctuation in the low energy limit. 

In order to analyze the superconductivity, the singlet 
pairing interaction is described by the following equation: 

Vsiq) = lsxs{q)s - ^CxMC +^iS + C). 

Then, the linearized Eliashberg’s equation is given as 

'X(t>hh{k - q)Gi3isiq - k), 

where (j)(ci) is the gap function. In the present study, the 
dominant form of the gap function has the s± form with 
some nodes on the electron Fermi surface for large phos¬ 
phorous content, consistent with previous studie o^^i^^~ — . 
The eigenvalue A monotonically increases with lowering 
the temperature T, and reaches unity at T = Tc. We cal¬ 
culate this quantity for a fixed temperature T > Tc, so 
that A(T) is positively correlated with, and thus measure 
of the Tc- It is known that the spin fluctuation around 
the wave vector (tt, 0)/(0, tt) contributes to the s± su¬ 
perconductivity, so in the actual calculation of r^(w), 
the summation over the wave vectors is restricted to the 
vicinity of (tt, 0)/(0, tt). 

III. RESULTS 
A. Bond angle dependence 

In Fi^ we show, as functions of the Fe-As-Fe bond 
angle a, the calculated superconducting eigenvalue A as 


well as the imaginary part of the low energy derivative 
of the intraorbital spin susceptibilities and 
A is calculated for s-wave symmetry, and the resulting 
gap function has the s± wave form. In ref, [l^ . some of 
the present authors have shown that the variation of A 
against the bond angle a exhibits a double dome feature 
as in Figl^l which is consistent with the experimental 
observation s . There, the origin of the double dome 
structure has been explained as follows: superconductiv¬ 
ity is suppressed when the bond angle is increased from 
a ~ 110 because the hole Fermi surface around (7r,7r) 
originating from the dxy orbital vanishes. 

However, the spin-fluctuation and A are re-enhanced 
in the larger bond-angle region (intermediate phospho¬ 
rous content regime) because the density of states of the 
inner hole Fermi surface originating from the d^z/yz or¬ 
bital increases. As the bond angle is further increased 
(near the phosphide end), the orbital character of the 
hole Fermi surfaces have the dxz and dyz character, 
where X-Y axes are rotated by 45 degrees from x-y. The 
electron Fermi surfaces around (tt, 0)/(0, tt) always have 
the dxz and dyz character, so that the matching of the or¬ 
bital character between the electron and the hole Fermi 
surfaces is not good for too large bond angles, thereby 
suppressing again the superconductivity. 

At first glance, the superconducting eigenvalue and the 
strength of the low energy spin fluctuation appear to be 
correlated. However, if we look more closely, we find 
that this is not necessarily the case around the first dome. 
Namely, (i) while A barely changes around the bond angle 
where it is maximized, i.e., 109 < a < 111, the strength 
of the spin fluctuation in that region strongly varies. At 
a = 109, on the left side of the dome, the dxz/yz spin 
fluctuation is large, but it is strongly suppressed at a = 
111, on the right side. Similarly, the dxy spin fluctuation 
is largely suppressed at a = 111 compared to a = 110 ^ 
109. (ii) On the left and right sides of the first dome, 
a = 106 and A = 112 have almost the same value of A, but 
the low energy spin fluctuation, especially that for the 
dxy orbital, is much larger for a = 106. (hi) At a = 109, 
the dxz/yz spin fluctuation exhibits a sharp enhancement, 
but this does not seem to affect the superconductivity. 
The low energy spin fluctuation is mainly governed by 
the details of the Fermi surface nesting condition, so this 
result indicates that those details and superconductivity 
are not intimately correlated. 

On the other hand, A around the second dome seems to 
be correlated with the strength of the dxz/yz low energy 
spin fluctuation. In this large bond angle regime, the 
dxy spin fluctuation is always suppressed because of the 
absence of the dxy Fermi surface around (tt, 


B. Temperature dependence of the low energy spin 
fluctuation 

In order to further investigate the correlation between 
superconductivity and the low energy spin fluctuation, 
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FIG. 2. Solid red:the eigenvalue of the linearized Eliashberg’s 
equation ; blue dotted and green dashed :the low-energy spin- 
fluctuation and ^ respectively; all against Fe-As-Fe 

bond angle a. We call the local maximum of A around a = IIQ 
and a = 118 the first and second domes, respectively. 


we now study the temperature dependence of 7 '^*' and 
ryxz/yz yg^j-jQ^g bond auglos. In FiglH] and FiglU we 
show the results for the bond angle within the first and 
the second domes, respectively. The development of both 
ryxzjyz ^xy jg yreak even at the bond angle of a = 
111 , where the superconducting eigenvalue is close to its 
maximum value. This further confirms our viewpoint (i) 
mentioned in the previous section. 

As for the second dome, increases upon lower¬ 

ing the temperature only at a = 118, where A is locally 
maximized . In short, the low energy spin fluctuation and 
superconductivity is not intimately correlated around the 
first dome, while the superconductivity is correlated with 
the low energy the spin fluctuation in the second 

dome. 


C. Energy dependence of the spin flnctuation 


So far we have seen that superconductivity is not 
tightly correlated with the low energy spin fluctuation 
around the first dome. The reason for this is expected 
to be because the finite energy spin fluctuation also con¬ 
tributes to the superconductivity. If the strength of the 
low energy and finite energy spin fluctuations were cor¬ 
related, so would superconductivity and the low energy 
spin fluctuation. Hence, it is expected that the low and 
finite energy spin fluctuations exhibit different behavior 
upon varying the bond angle. To see this, we go back to 
the issue raised as (ii) in section ITlIAl and compare r(a;) 
between the bond angles a = 112 and ol = 106 by taking 
their ratio as in Fig IS] It is found that at some finite 
energies, F(a;) at a = 112 is larger than that at a = 106, 
namely. 



FIG. 3. The measure for the low-energy spin fluctuation 
r^xz/yz r^xy (^bottom) aguiust tcmperature for vari¬ 

ous bond angles from 108 to 112 around the first dome. 
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FIG. 4. Plot similar to Fig[3] for the bond angles 116, 118, 
and 120 around the second dome. 
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FIG. 5. The ratio of (top) and (bottom) 

between bond angle 112 and 106 plotted against uj. Solid red 
(dotted blue) lines indicate regions where the spin fluctuation 
for a = 106 dominates over that for a = 112. 


A(106) ~ A(112) 

7a=106 > 7a=112 
^0=106 (^finite) (^finite) 

This confirms our expectation that the absence of the 
tight correlation between superconductivity and the low 
energy spin fluctuation in the first dome is due to the con¬ 
tribution from finite energy spin fluctuation, which itself 
is not correlated, or in some cases even anti correlated, 
with the low energy spin fluctuation. 

In the second dome, on the other hand, the magnitude 
relation of r(a;) among the bond angles a = 116, 118, 120 
depends on oj as seen in Fig[Hl The only clear correlation 
between spin fluctuation and superconductivity is seen 
at very low frequencies. Therefore, for the second dome, 
the low energy spin fluctuation seems to dominate the 
variation of the superconducting eigenvalue A. 

IV. DISCUSSION 

As described in the Introduction, it was found in the 
NMR experiments^ii^ that the correlation between Tc 
and the development of 1 [TiT at low temperatures is ma¬ 
terial dependent. In ref.[ 8 | in particular, it has recently 
been found that the enhancement of in the interme¬ 
diate regime of the As content in LaFeAs^Pi-xOi-yFy 
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FIG. 6. Plots similar to Fig[5] The ratio of between 

a = 116 and a = 118 (top), and a = 120 and a = 118 
(bottom) 


is correlated with the development of l/TiT at low tem¬ 
peratures. On the other hand, the development of 1/TiT 
tends to be saturated upon lowering the temperature in 
a material (Y^Lai-zFeAsOy) reaching Tc = 50K, where 
La is partially substituted by Y in LaFeAsOy. 

These experimental results can be understood in terms 
of the calculation results obtained in the present study 
by notifying that (i) 7 ^^ and mainly contribute to 

l/TiT, because xy^ yz, xz orbitals are the main origin 
of the spin fluctuations, and (ii) increasing the phospho¬ 
rous content corresponds to enlarging the bond angle a, 
while substituting La by Y corresponds to reducing the 
bond angle. LaFeAsOi-xF^ has the bond angle of about 
a = 113 ^ 114, so Y^Lai-zFeAsOy is expected to be 
sitting in the first dome, where the dxy spin fluctuation 
plays an important role, and LaFeAsjjPi-xOi-yFy is in 
the second dome, where only the d^ziyz spin fluctuation 
is important. As was mentioned in the previous section, 
7 tends to saturate or even decreases upon lowering the 
temperature in the right side of the first dome (bond an¬ 
gles a = 111 ^ 112 ) even though the superconducting 
A values are substantially larger than those in the sec¬ 
ond dome. On the other hand, the local enhancement 
of A is correlated with the development of at low 

temperatures in the second dome. 

To see more directly the correspondence between the 
calculation and the NMR experiment, we plot in Figl3 
a summation 27 “^/^^ -|- 7 ^^/^^, which should give the 
dominant contribution to NMR l/TiT, renormalized by 
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FIG. 7. + 7 “^^/^^ normalized by its value at T = 

O.OSeV, plotted against T for a = 111.5 and 118. 


V. CONCLUSION 

In the present work, using the five orbital model of the 
1111 iron-based superconductor, we have studied the cor¬ 
relation between the spin fluctuation and the supercon¬ 
ducting transition temperature (eigenvalue of the Eliash- 
berg eq.), which exhibits a double dome feature upon 
varying the Fe-As-Fe bond angle. Around the first dome 
with higher Tc, the low energy spin fluctuation and Tc 
are not tightly correlated because the finite energy spin 
fluctuation also contributes to superconductivity. In fact, 
even near the Tc maximum, the measure of the low energy 
spin fluctuation, proportional to 1 /TiT, tends to saturate 
or even decrease upon lowering the temperature. On the 
other hand, the magnitude of the low-energy spin fluctua¬ 
tion originating from the d^z/yz orbital is correlated with 
Tc in the second dome. These calculation results are con¬ 
sistent with recent NMR study^, and therefore strongly 
suggest that the pairing in the iron-based superconduc¬ 
tors is predominantly caused by the multi-orbital spin 
fluctuation. 


its value at T = O.OSeV for angles a = 111.5 (in the 
right side of the first dome) and a = 118 (at which A 
is locally maximized in the second dome). Indeed, this 
quantity starts to decrease for T < O.OleV for a = 111.5, 
while it continues to develop down to low temperature for 
a = 118, in nice correspondence with the experiment in 
ref. Q . These results strongly indicate that the spin fluc¬ 
tuations originating from multiple orbitals are the origin 
of the superconductivity in these materials. 
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